function [linVals MorletVals] = evalMorlet(LB, UB, N)
  dx = (UB-LB)/(N-1);
  for k=1:N,
      cv = LB+k*dx;
      linVals(k) = cv;
      MorletVals(k) = morletVal(linVals(k));
%      MorletVals(k) = exp(-(cv*cv)/2)*cos(5*cv);
  end
  
  % LB = -4; UB = 4; N=30000; tic; [lv mv] = evalMorlet(LB, UB, N); tm=toc; plot(lv, mv); 
  % tm = 1.58sec